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Liischer's "admissibility" condition on the gauge field space plays an essential role in constructing 
lattice gauge theories which has exact chiral symmetries. We apply the gauge action proposed by 
Liischer with the domain-wall fermion action to the numerical simulation of the massive Schwinger 
model. We find this action can generate configurations in each topological sector separately without 
any topology changes. By developing a new method to sum over different topological sectors, we 
calculate meson masses in nonzero-6 1 vacuum. 
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I. INTRODUCTION 

There are various problems in gauge theories where the chiral symmetry plays a crucial role. Although the lattice 
gauge theories provide a method for nonperturbative computation on these problems, they have not given satisfactory 
answers since the conventional fermion action suffers from either the species doubling 0, 0] or the lack of chiral 
symmetry, which makes the study of the chiral behavior difficult. For this reason, extensive studies have been made 
to improve the fermion action. It was found that lattice Dirac operators satisfying the Ginsparg- Wilson relation Q, 
which is for example realized by the Neuberger's overlap Dirac operator 0,0, have exact chiral symmetries 0. Such 
actions are expected to give a fundamental improvement in the study of K meson physics, finite temperature QCD, 
or even chiral gauge theories, despite their complicated forms. 

At the classical level, the Ginsparg- Wilson relation is sufficient to solve the problem of chirality. However, at the 
quantum level, topological properties of the gauge fields in the continuum space should also be kept on the lattice in 
order to reproduce the correct chiral anomalies. Liischer found that one can construct lattice gauge theories without 
breaking topological structures by restricting the link variables to satisfy the "admissibility condition" || ; 

||1 - Pfiu{x))\\ <e for all x,[i,v. (1) 



where 



P„ v {x) = U^(x)U v (x + fflUtix + va)Ul{x), (2) 



and e is a fixed positive number. In order for the gauge fields to satisfy this condition automatically, he proposed the 
following action as an example, 



a (1 - RoP M „(x)) .. , ... 

^H'.S. '-IIi-pZwIIA 1 " 11 -^" 1 ^. (3) 

oo otherwise 



The admissibility condition makes gauge fields smooth and unphysical configurations such as vortices are sup- 
pressed. In fact, the space of admissible fields is separated into disconnected subspaces labeled by some integers which 
correspond to topological charges in the continuum theory 0,^3- Therefore, one can precisely treat topological effects 
such as U(l)problem, 8 vacuum and so on. Moreover, Liischer's action is also differentiable and gauge invariant and 
has a good continuum limit as Wilson's action. 

Although there are other proposals for improvements of the gauge action [Til Il2l Il3l IT3. ITEl Il6l \vi\ , they do not 
keep any topological structures. In this sense Liischer's action combined with Ginsparg- Wilson fermion action would 
be the best choice to investigate chiral symmetries on the lattice. Moreover, admissibility condition is indispensable 
to construct chiral gauge theories since the gauge symmetry is never realized without exact chiral symmetries Q . 
Thus it would be important to examine how admissibility works and how topological structures are realized on the 
lattice in numerical simulations. 

In this paper, we apply Liischer's gauge action to the numerical studies of the two-flavor massive Schwinger model 
|19| on the lattice using the domain wall fermion action pfl l2f| in order to examine how the admissibility works by 
probing the topological properties of the lattice theories. The massive Schwinger model is a good test ground for 
several reasons; it has been well examined analytically in continuum space even in strong coupling limit, its vacuum 
has non trivial topological structures due to the chiral anomaly, and it also shares many other interesting properties 



2 



with QCD such as U(l) problem and confinement. There already exist extensive studies of to pologic al structures of 
the massive Schwinger model have on the lattice with Wilson's gauge action in the literature 

where our work provides an alternative lattice approach with Liischer's action. Although it will be interesting subject 
to make a detailed comparison of the results with our method with previous calculations, we will leave it for the future 
publication and focus on the feasibility study of our method and the new results on 9 vacuum in the present work. 

We found that Liischer's gauge action can generate configurations in each topological sector separately. We develop 
a new method to evaluate the observables in nonzero-0 vacuum by summing over those in different sectors with 
correct wei ghts . (Our strategy is quite different from that of sampling other sectors by enhancing topology changes 
[29l l3fl I31U32II.) We applied our method to the meson correlators and observe the dependence of the isotriplet 
meson mass. We reproduce the well-known continuum results; the isotriplet meson mass scaling as a function of the 
fermion mass and the their 9 dependence, and the fact that the isosinglet meson acquires a heavier mass than the 
isotriplet meson due to anomaly (the so-called U(l) problem). 

This paper is organized as follows. In section [HJ we summarize main results of the continuum massive Schwinger 
model. In section lnTl we discuss details of our simulation. In section lTVl we present the results and compare them with 
continuum theory. In Appendix, we compare Liischer's action and Wilson's action in the quenched approximation, 
and show the impact of the admissibility condition on the topological and chiral properties. 



II. REVIEW OF THE MASSIVE SCHWINGER MODEL 



A. Continuum theory 



We consider the two flavor massive Schwinger model |33u34L l35| with degenerate masses. The continuum action in 
Euclidean space is defined as 

S = Sq + Sf, 
Sg = J d 2 x^F^(x)F^(x), 

S F = f <fxJ2^){P+m)^(x), (4) 

i—l 

where 

2 

F^(x) = d„A v {x) - d v A„(x), p=Y,^(dn + iA„), 

V-f ? nV 7 2 =fr;V 7 3 = -7V, (5) 



is the gauge field and tp is the two-spinor fermion field. We take g and m to be positive without losing generality. 
If we take the space-time to be a torus T 2 , the space of gauge fields is separated into topological sectors each of 
which is labeled by an integer 



1 



4tt 



N=— d 2 xe^F^, (6) 



T 2 



where we take the sign convention of the antisymmetric tensor as 612 = 1. Then this theory has vacuums dependent 
on phase 9. Full path integrals are defined by a summation of integrals in each topological sector; 



+00 



Z full (g,m,6) = e lN9 ZN(g,m) 

N=~oo 

Z N {g,m) = J DA% DrfDipexpi-Sa - S F ), (7) 

where A^ denote gauge fields in the sector with topological charge N. Using Eq.©, we can rewrite Zf u u(g, m, 9) as 
follows 

Z fuU (g,m,9) = J DA^DiPeMSG - S F - S e ), (8) 
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where 



S e = -i I d 2 x—e^F^. (9) 



It is well-known that this model is equivalent to the two-component scalar theory 34, 35, 3£ 

,2 



S = d 



-j2 , 



(10) 



— 2cmgcos I v 2n(j)+ (x) — — J cos(v 2tt(/)-(x)) 
where /io — g\ — and c is a numerical constant. 

V 7T 

For m <C /io an d 0^0, perturbative calculations of 0(m) show that light scalar tfi_ has the mass 

m- = V2^{2cm^l 12 cos ^) 2/3 , (11) 

Li 

and heavy scalar <p + has the mass 

m+ = A*o + {0{m) corrections). (12) 

Now, we discuss the chiral behaviors of the two-flavor Schwinger model in the chiral limit m — > 0. The action 
has U(2)a — SU(2)a X U(1)a chiral symmetry in this limit. U(1)a symmetry is broken by the anomaly, which 
manifests itself in the vacuum with nontrivial topological structures. We define the isotriplet meson operator t:q = 
^173^1 — V'273'02 and the isosinglet meson operator rj = ^173^1 + ^273^2 by the fcrmion bilincars. In the bosonization 
picture, it is shown that ttq propagation corresponds to that of light scalar 0_ and rj propagation corresponds to that 
of heavy scalar <fi + . In the massless limit, Eqs. <|11[) and (|12fl show that ttq becomes massless while r\ remains massive 
in accordance with the U(l) problem in two-dimensional QED. 

B. Lattice theory 

Let us consider the lattice regularization of the massive Schwinger model. We take the lattice size is L x L x L3 
with lattice spacing a = 1, where L3 is the length of the third direction for the domain wall fermions [2fT. I2H] . Action 
is defined as follows 

S = pS G + S F , (13) 
(l-ReP^x)) 



; IT 1" y ,' '— , if admissible 

1 - (1 



S G = { ^l-(l-ReJV(a:))/e , (14) 

00 otherwise 
2 

Sf = J2T,T,[^ d dw(x,s;x', s 'WAx') 

x,x' s,s' i—1 

+<f>l*(x)D AP (x,s;x',s')<Pl,(x')], (15) 



where 



! 2 

D DW (x,s;x',s') = ^^{(l + l^U^S.. 



+ (1 - Jn)Ul(x - A)^-A,x'^s,s' } 

+(M-3)5 x , x ,5 StS , 

s' u x,x' 

+(m - l)P + 6 SiL3 5 s > A 5 XtX , + (m - l)P-S Si iS s > tL3 6 x . x >, 
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1 2 

D AP (x,s;x' ,s') = - {(1 + Jn)Un(x)S x+ n >x t6 a , a i 

+ (1 - 7 Al )[// t (x - fx)S x -p, tX 'S SiS >} 
+(M - 3)5 x , x ,S SiS < 

+ P+3s+l,s'Sx,x' + P-fis-l,s'$x,x' 

-2P + 5 St L 3 5 s > t i5 XjX i - 2P-5 sA 5 sl , L3 5x,x'- (16) 

(i = 1/g 2 , M is a constant satisfying < M < 1, denotes summation over all plaquettes and P± arc the chiral 

p 

projection operators; 

P ± - ^ (17) 

m is the fermion mass. </> l 's are Pauli-Villars regulators which cancel the bulk contribution. 

Since it is not possible to change the topological charge by a local updation under the admissibility condition, our 
lattice theory with Liischer's gauge action has a topological invariant 

N = -^^lnP 12 (x). (18) 

X 

This charge corresponds to Eq.@ and gauge field configurations are classified into topological sectors . Each sector 
characterized by N has the classical gauge configuration U^(x,y) minimizing the action, which is given as 

U? [N \x,y) = exp j h^LL - ZZ^SxX-ix } , 



1-KWx 


2nNi 


L 


L 




2-nNi 


L 


+ V 



uf [N \x,y) = expj^p + ^yj, (19) 

up to gauge transformations, where v\ and are the parameters which determine the values of Wilson lines in x 
and y directions. v\p, can take any values in the region < v\ 2 < 1 • This configuration gives constant background 
electric fields over the torus. 



III. LATTICE SIMULATIONS 



A. Observables in each sector 



Simulation is carried out by Hybrid Monte Carlo method with Liischer's gauge action Ea. ll4|) . The matrix inversions 
are calculated by the conjugate gradient algorithm. 

We take 16 x 16 x 6 lattice. The simulation is carried out at (3 = 1/g 2 = 0.5 and M = 0.9. The parameter for 
the admissibility condition is chosen as e = 1.0. At this value of e, we find that initial topological charge is not 
changed through the simulation. ( See Fig. n an d Fig. [3 ) For the fermion mass, we choose m = 0.1,0.2,0.3,0.4. 
50 molecular dynamics steps with stepsize At — 0.02 are performed in one trajectory of the Hybrid Monte Carlo 
algorithm. Configurations are updated per 10 trajectories. We generate 500 configurations for each topological sector 
by taking the classical configuration in Ea. (|19fl as the initial configuration. From the set of configurations in each 
sector with topological charge N, we measure the isotriplet meson propagator 

O^) =]>>(*, y)7r(0,0))£; m , (20) 
y 

and the isosinglet meson propagator 

y 

where ()? m denotes the expectation value in the TV sector. 
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FIG. 1: The comparison of the Monte Carlo evolution of the topological charge with Wilson's gauge action and Liischer's gauge 
action for the same lattice spacings determined from the string tension. Initial topological charges are zero. Left: Wilson's 
gauge action with (5 = 2.0, Right: Liischer's gauge action with (5 — 0.5. The topological charge changes for Wilson's gauge 
action, while it does not for Liischer's action. 
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FIG. 2: The Monte Carlo evolution of the topological charge with Luscher's action for (5 — 0.5 is shown. The initial topological 
charge is two. 



B. A new method of summing over different topological sectors 



The Hybrid Monte Carlo simulation is performed by small changes of link variables. Thus choosing the configuration 
given by Eq. i|19|) as initial condition, we can generate configurations without changing the topological charge for any 
value of the coupling constant. 

Now, we discuss full path integrals on ^-vacuum. Suppose that we measure the expectation value of an operator O, 



/n\' ul1 — £-/N=—oo J g y y / 00 \ 



where U„ denote link variables in the sector with N and 



Z N (P,m) = J DU™ DipDipe~ p G ~ F (23) 

is the lattice counter part of the Zn in Eq.JJJ). In terms of the expectation values in each topological sector, (O)^^ 
can be rewritten as, 

(0) fuu = E + N =-ooe me (0)l m R N (P,m) 
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where 



and 



(G) N 



{3,m 



J DU^ DipDip Oe-t 3SG ~ s F 
Z N (f3,m) 

Z N (j3,m) 



R N (/3, to) 



Z ((3,m) 



We call R (P,m) the reweighting factor. Note that Z N {(3,m) satisfies the following differential equation; 



dZ N (J3,m) 
d/3 



/Z N {(3,m) 



~(Sg) fam- 



ily integrating over (3 again, Z^{f3,m) is expressed as, 

Z N (/3,m) = Z N (oo,m)exp (^J d(3' {S G )p, ^ . 



(25) 



(26) 



(27) 



(28) 



Then, the reweighting factor R N (/3, to) is expressed as, 



R N (/3,m) 



i?/v(oo, to) 
Z (oo,m] 



■ exp 



((S G ' N 



■>G)f3> 



exp(— (3S, 
x exp 



N 



J d^ 2 det(D^ w ) 2 /det(D^ P ) 



G mmJ f dv x dv 2 det(D° DW ) 2 / det(D° AP ) 2 







(29) 



where S^ mm 



is minimum of the gauge action in sector with TV given by constant background fields Eq. (|19(l and D^ w 



and D^p are Dirac operators given by this background. Note that the integrand vanishes rapidly as j3' — > oo, so the 
integral over (3' converges. 

We can evaluate the full path integrals on 0-vacuum by obtaining (0)p m and R N (f3, m) in each sector. It should be 
noted that this method is only possible with a gauge action with the admissibility condition in which the topological 
charge is strictly conserved. In our approach we take the property that Liischer's gauge action allows no topology 
change at all as an advantage and treat the sum over the topologies in a controlled fashion. This is in contrast to the 
conventional gauge actions, with which the topology change is suppressed but not completely prohibited so that one 
has to tackle the problem of enhancing topology changes. 

A related but somewhat different approach was proposed by Diirr j3]j where one makes a quenched calculation 
and give the whole fermion determinant as the reweighting factor. He also proposed an approximation in which one 
replaces the determinant for the given configuration by the determinant of a common representative configuration for 
the given sector, which reduces the enormous computational effort . 

Of course, computing the reweighting factors for the sum over topologies requires extra works. Whether this 
program works must be examined in practical simulations. In the following subsections we show that our new method 
is valid and full path integrals can be evaluated with controlled statistical and systematic errors. 



C. Calculation of R N (/3,m) 

Let us discuss how to evaluate R N (f3,m). The classical minima of the gauge action S G min are evaluated easily. In 
Fig. |3 we can see Sq min is numerically proportional to \N\ 2 . 

The fermion determinant det D 2 on classical background in the sector with ./V is numerically calculated using the 
Householder method and the QL method [3^|. The integral over the moduli v\p is approximated by the weighted 
sum over the discrete set of points uniformly distributed in the whole integration region as Fig. The numbers of 
points for the weighted sums are 5 x 5 for both det(D° DW ) 2 / det(L>^ P ) 2 , and for det(D% w ) 2 / det(D% P ) 2 with N ^ 0. 
The value of 



= / dvidu 2 det(D% w ) 2 / det(D% P 
f d Vl dv 2 det{D% w ) 2 1 det{D% 
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is plotted in Fig[5] It decreases as \N\ increases, due to the contribution of small eigenvalues proportional to the 
fermion mass, which emerge in the nontrivial topological sectors since Atiyah-Singer index theorem is realized on the 
lattice in 273—5-00 limit [7j. 

In order to obtain the exponential factor in Ea. l29fl . we need to evaluate the integral of the following quantity 

Ssubtrifl'i m ) = {&G ~ &G min)p',m ~~ {^g)% ,m- (31) 

Since S^ ubtr {(3' ,m) decreases rapidly as (3' — > oo, the integral of S^ ubtr over f3' is well approximated by a weighted 
sum over the discrete set of points for S^ ubtr (f3' , m) at (3' = 0.5, 1.0, 1.5, 2.0. For each (3', we evaluate S^ ubtr (f3' , m) by 
sampling more than 5000 configurations. 

Total reweighting factor R N ((3,m) at p = 0.5 and m = 0.2 is plotted in Fig. [S] It is shown that higher topological 
sectors are indeed suppressed by the reweighting factor. 

Finally, combining the correlators and the reweighting factors, we obtain the total expectation values on nonzcro-0 
vacuum as 

4 

5>( a M,M0 J 0)) /trfI = e im J2«^vMO,O))^ m R N (0,m), (32) 



V JV=-4 y 

and 

4 

J2(^,yM0,0))full= E elN6 E^^' 0))l m R N (f3,m), (33) 

y JV=-4 y 

up to a constant normalization factor. Here we have ignored \N\ > 4 sectors since they only give contributions less 
than 1.2 % of zero sector. Then we can get pion mass and r\ meson mass including full nonpcrturbative effects and 9 
dependence. In this calculation, the propagators are fitted by minimizing the x 2 an d the total statistical errors are 
estimated by summing those in individual sections in quadrature. 
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FIG. 3: Minimum action in each topological sector is plotted as a function of the topological charge. Filled squares are the 
data and the dotted line is the fit with a quadratic function. 
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FIG. 4: Three dimensional plot of the v dependence of the fermion determinant det(Dp W ) 2 / det(D^ P ) 2 . Left: N = case, 
Right: N = 1 case. 
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FIG. 5: N | dependence of Det N for the fermion mass m 
= 0.1, 0.2 are shown. 
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FIG. 6: The total reweighting factor #^(0.5,0.2) is plot- 
ted as a function of the topological charge. The factor 
falls off rapidly as the topological charge increases. 



D. Systematic errors 



In this subsection we discuss possible systematic errors. These error estimations show that our simulation is 
reasonable and the results are reliable. 

Let us now study the lattice spacing dependence. We measure a dimensionless quantity A = m^/(m 4 a) at (3 = 
0.5, 1.0, 1.5,2.0 in zero topological sector, where m w denotes pion mass and a denotes the string tension. As Fig. 
indicates, the results at (3 = 0.5 show no large lattice spacing dependence which suggests that the discretization error 
is under control. 

Next we discuss finite size effects for the space-time size L and and for the extra dimension size L 3 . We measure 
the pion mass on the lattices of size L 2 x L 3 = 8 2 x 6, 10 2 x 6, 16 2 x 6, 20 2 x 6, 16 2 x 2, 16 2 x 4 and 16 2 x 10 in the zero 
sector. Fig. [5] shows L dependence and Fig. [5] shows L3 dependence. We find the meson mass is stable for L larger 
than 16 and for L3 larger than 6 so that the finite size error is also under control with our choice of the lattice size 
16 2 x 6. The discretization error and finite size errors from the nonzero topological sector is similarly under control. 

We now study the error in the integration over the moduli v\p. In order to estimate the systematic error we 
also evaluate the integral by the weighted sum of 10 x 10 points. We find that the change is very tiny ( relative 
change ~ 10 -8 ) and is negligible compared to other systematic errors, as is expected from the mild v dependence 
of det(D° DW ) 2 / det(D° Ap ) 2 in Fig. H Fig. H also shows that det(D% w ) 2 / det(D% p ) 2 with N ^ has almost no v 
dependence. In fact this remarkable flat dependence is also seen in the continuum counterpart analytically |39|. We 
therefore conclude that the error in the weighted sum is even more negligible for the nonzero topological sector. 

Since the integral of S^ ubtr over 0' is approximated by the trapezoidal rule using the data for the discrete set of /?' 
points, the error in this approximation should be estimated. For this purpose we evaluate the integral of <S'^ btr in an 
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By trapezoidal rule 


By fit 


.R (0.5, 0.2) 


1.0 


1.0 


pi/n n o\ 
JrC (U.o, U.z j 


U.OO ( \oK) ) 


U.oy^zz ) 


R 2 (0.5, 0.2) 


0.442(39) 


0.45(16) 


i? 3 (0.5,0.2) 


0.201(25) 


0.32(18) 


i? 4 (0.5,0.2) 


0.0636(91) 


0.072(46) 


(at = 0) 


0.647(07) 


0.650(34) 



TABLE I: The reweighting factor in each sector from two different methods of evaluating the integral of S^ ubtr ; The trapezoidal 
rule and the integral of the polynomial fit. The resulting pion masses are also given. 



alternative way, in which we fit the discrete set of data with the function of the form 

S^ubtrW' m ) — ~ 72 + ~jj 73 ' (34) 

and compute the integral of S^ ubtr over j3' analytically. Table [I] shows the difference of the two ways of evaluation. 
The pion mass are consistent with each other. Thus we find that the approximation for the integral of Sf ubtr does not 
give large systematic errors in the meson mass. We now study the truncation error in the sum over topological sectors. 
As we discussed before, we neglect | | > 4 sectors since these contribution are suppressed by large value of action 
and fermion zero modes. Fig. I1UI shows the pion mass at 9 = 0.3n measured for a variety of the highest topological 
charge N max . Therefore the truncation error in the sum over topological sectors are negligible in comparison with 
the statistical errors. 
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FIG. 7: Lattice spacing dependence of the dimensionless quantity A = ml\/(m 4 a). Horizontal axis is the string tension in 
lattice unit. 



10 



data 
0.609 



data 
0.609 



0.8 



0.8 



0.6 



0.6 



0.4 



0.4 



0.2 



0.2 



I 1 1 1 1 1 

5 10 15 20 25 

lattice size L 

FIG. 8: The dependence of the pion mass on the lattice 
size L in the space-time direction. The fermion mass is 
m = 0.2. Filled symbols are the data and the dotted line 
shows the fitted mass for L = 16. The pion mass shows 
no volume dependence for L > 16. 
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FIG. 9: The dependence of the pion mass on the lattice 
size I/3 in the third direction. The fermion mass is m = 
0.2. Open symbols are the data and the dotted line shows 
the fitted mass for Lz = 6. The pion mass shows no 
volume dependence for Lz > 6. 
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FIG. 10: The pion mass from the sum over \N\ < |7V„ 
mass for 7V max > 4 is observed. 



sector contributions for m = 0.2 and 8 = O.Stt. No change in the 



IV. MESON MASSES 



A. Pion mass and 8 dependence 

Fig. ^2 shows pion propagators in each topological sector and Fig. ^| shows full propagators at various 8. We 
measure the pion mass by fitting these data to the hyperbolic cosine function. The fit range is x = [5, 8] for which we 
find good plateau in the effective mass plot as Fig. 1131 shows. In fitting x 2 /dof is also a small value (x 2 /dof < 0.1). 

Figll4l shows pion mass at 8 = as a function of fermion mass m. We ignore the m dependence of S^ ubtr (f3' ,m) 
and use m = 0.2 result for all m. We fit the results to the following function suggested by the continuum theory with 
possible additional constant term b from the residual mass of pion, 

m„(m) = am 2/3 + b. (35) 

Figll4l shows that Eq. (|35[1 fits the data very well (x 2 /dof = 0.39) so that the fermion mass dependence is consistent 
with the continuum theory. The residual mass of the pion measured in the chiral limit is also tiny as b — — 0.057±0.060, 
which shows that the violation of the chiral symmetry is very small. 

In Fig. E| we present the 8 dependence of the pion mass at (3 — 0.5 and m = 0.2. As a remarkable feature, the 
result is in perfect agreement with that in the continuum theory in 0/(2-71") < 0.5 region. A good control of the 8 
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dependence shows that our method for summing over different topological sectors with Liischer's gauge action indeed 
works numerically. 

At large 9, statistical errors increase, which are due to cancellations of propagators among different topological 
sectors. In the calculation, we approximate the integral of Sf ubtr (f3' ,m) by the trapezoidal rule for the discrete set 
of j3' points, but this does not seem to be the reason for the large fluctuation in the 0/{2it) > 0.5 region. The main 
nonperturbative contribution comes from Det N and S^ ubtr ((5' , m) gives only perturbative effects of order (3'~ 2 . In 
fact, even ignoring the integral (we set S^ ubtr (f3' , m) — for all /?') we can get similar results as Fig. 1101 

We suspect that this large fluctuation is an example of the well-known phase problem. Simply increasing the 
statistics might not improve the situation. 

Of course in application to QCD, it will be important to evaluate S^ ubtr (f3' , m) and other observables more precisely. 
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FIG. 11: The pion propagator in each sector for m = 0.2. 



0.01 



e = 0.0 i — i- 

9=0.1 7i h-x- 

e =o.2k >•--•*• 

6 =0.371 i ED 



□ 
i 

10 



12 



14 




FIG. 12: The full pion propagators with m = 0.2 for FIG. 13: The effective mass plot of the pion for m = 0.2 
various 6 are plotted. and 9 = 0. The dotted line shows the result of the fit. 
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FIG. 14: The fermion mass dependence of the pion mass for 8 — 0. The crosses are the lattice data and the dotted line is the 
result of the fit with the function in Eg. 13511 . The chiral behavior is consistent with that of continuum theory. 
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FIG. 16: 8 dependence of the pion mass obtained by ignoring S 3ubtr (f3' , m) in the reweighting factors 



B. 77 meson correlator and U(l) problem 

As the final subject, we would like to present the result of our exploratory measurement of the 77 meson mass in 
order to study the topological structure. The 77 propagator consists of two parts; 



>/'/; : -1(1 >■[ ^3^73-^ 



Mi>- { 73^ 



tr 



73 -p )} 



(36) 



where the first term is the same as flavor non-singlet 7r propagator and the second term gives the "hair-pin" or 
disconnected contribution to the flavor singlet operator. Because the number of physical space-time points is only 
16 x 16, we compute the "hair-pin" contribution by brute force, namely by solving the fermion propagator for all 
points without relying on the noise method^O] or Kuramashi method |4lj . 

Fig. El shows the contribution of the second term in each sector, whereas Fig. ^5] shows the full (symmetrized) 
77 propagator at m = 0.2 and 8 — 0. We also present effective mass plot in FigllSI We find that the fall of 77 
propagator is steeper than that of ir which gives qualitatively consistent results with U(l) problem, although it suffers 
from large statistical fluctuations making quantitative studies difficult. Further quantitative studies require some new 
ideas which efficiently reduce the statistical error and controls the systematic errors due to the summation over large 
number of topological sectors. 
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FIG. 17: The propagator of 77 in each sector at m = 0.2. 



FIG. 18: The effective mass plot of the full r\ propagator 
for m = 0.2, 6 = 0. 
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FIG. 19: The full propagator r\ at m = 0.2 and 6 — ( filled squares ). The pion propagator is also plotted for comparison ( 
open circles ). The propagators are normalized by the value at x = 1. 



V. SUMMARY AND DISCUSSION 



In this paper, we elucidate the role of the admissibility condition on the topological and chiral properties in lattice 
gauge theories by applying Liischer's action together with domain wall fermions to a numerical simulation of the 
massive Schwinger model. To investigate the ^-dependence of the correlators, we have developed a method to sum 
over different topological sectors. We have found that Liischer's action is indeed applicable to Monte Carlo simulations 
and all the results are consistent with those in the continuum theory, confirming the validity of our method. 

We summarize the features of this action here again. (1) In Liischer's action, the gauge field strength is uniquely 
determined from the plaquette and the gauge action is a smooth function of the field strength. (2) The range of the 
action is not compact; 

< S G < oo. (37) 

This is the same situation as continuum theory We can treat the theory in terms of the field strength rather than 
plaquettes. According to these features, Liischer's gauge action has many advantages. 

1 . The use of this gauge action with the domain wall fermion action is valid even for the strong coupling regime since 
unphysical configurations are suppressed. (We find the suppression effect is especially remarkable in quenched 
approximation as discussed in Appendix.) 

2. We can treat the topological properties of the lattice theories precisely. This exact topological treatment is 
useful not only mathematically but also in a practical point of view. In conventional approach, there are two 
technical problems, i.e. violation of chirality at strong coupling and the slowing down of the topology change 
in unquenched simulation. For the former problem, the improved gauge actions which suppress the dislocations 
is proposed. However, in principle the suppression of the dislocations also suppresses the topology change so 
that the latter problem becomes even more difficult. Our method makes the improvement to the extreme and 
prohibits both the dislocation and the topology change completely, however by computing each topological 
sector and its reweighting factor we can reconcile the solutions to the the topology change problem and the 
dislocation problem at the same time. 

3. Once each topological sector can be computed separately, we can obtain a 9 dependence at once. 

4. Aside from the fact that we must simulate for each sector the typical simulation time needed for the trivial 
topological sector is no larger than that of using Wilson's plaquette action. For the nonzero topological charge 
sector, one can also increase the statistics at will very efficiently, in contrast to the conventional method where 
one can increase the statistics only by reaching the thermal equilibrium. In this sense, our method would have 
advantages in physical quantities for which the topological sectors with larger instanton numbers give larger 
contributions. 
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We expect that it would not be difficult to apply Liischer's type of gauge action also to QCD in four dimensions. In 
order to study the 9 vacuum, one should study how the reweighting factors can be computed. Also one should study 
the exact topological index on the lattice since the topological properties may be very complicated in 4-dimensional 
QCD on a torus. At least we can recommend the application to the calculation in the exact chiral limit since we 
usually need only zero sector results. We hope the understanding of the topological properties in lattice QCD will be 
improved by applying Liischer's admissibility condition. 

Finally we would like to point out that the method proposed in this paper to use Liischer's type gauge action and 
sum over topological sectors would also be essential for the numerical simulations of the chiral gauge theories on the 
lattice in the future. 
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Appendix 

In this appendix, we examine the validity of Liischer's action in the quenched approximation. In the strong coupling 
region, we compare Liischer's action with Wilson's action ignoring fermion loops. We set lattice size to be 32 x 32 x 5 
and measure pion mass in the zero sector. Gauge coupling f3 is chosen to give the same string tension a = 0.18 ; 
/3 — 1.0 for Liischer's action and j3 = 3.4 for Wilson's action. 

Fig l20l shows evolution of topological charge. In quenched approximation with Wilson action, fermion zero modes 
are all neglected so there is no suppression on topology changes. As a result, much unphysical configurations are 
generated. On the other hand, Liischer's action never allow topology changes. In Fig[^] the difference is clear. 
Liischer's action gives consistent pion mass and good chiral limit even in quenched approximation at very strong 
coupling. 

For a theoretically complete study of the quenched Schwinger model, we should take a sum over different topological 
sectors. However in our "quenched" study , unlike in the unquenched case, only the sector with zero topological charge 
is taken for the calculation with Liischer's gauge action, in order to illustrate the effect of the admissibility condition 
with the simplest example. It will be necessary to perform a theoretically complete quenched study by summing over 
different topologies and compare with analytic result in which it is predicted that the quenching effect does give a 
different fermion mass dependence from that in the unquenched theory |27| . 
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FIG. 20: The Monte Carlo evolutions of the topological charge in the quenched calculation with Wilson's gauge action and 
Luscher's gauge action for the gauge couplings having the same string tension. Left: Wilson's gauge action at f) = 3.4. Right: 
Luscher's gauge action at f3 = 1.0. Luscher's gauge action shows no topology change. 
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FIG. 21: The chiral behaviors of the pion mass in the quenched calculation with Wilson's gauge action and Luscher's gauge action 
for the gauge couplings having the same string tension. In this quenched study for illustrating the effect of the admissibility 
condition, only the sector with topological zero is taken for the calculation with Luscher's gauge action. Left: Wilson's gauge 
action at f3 = 3.4 . Right: Luscher's gauge action at /3 = 1.0 . Wilson's gauge action suffers from large fluctuation while 
Luscher's gauge action shows a good chiral behavior. Both of them are calculated by domain-wall fermions. 
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